Synthesis and structural characterization of 2Dioxane-2H 2 OCuCl2: metal-organic 
compound with Heisenberg antiferromagnetic S=l/2 chains 
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A novel organometallic compound 2Dioxane-2H20-CuCi2 has been synthesized and structurally 
characterized by X-ray crystallography. Magnetic susceptibility and zero-field inelastic neutron 
scattering have also been used to study its magnetic properties. It turns out that this material 
is a weakly coupled one-dimensional S=l/2 Heisenberg antiferromagnetic chain system with chain 
direction along the crystallographic c axis and the nearest-neighbor intra-chain exchange constant 
J=0.85(4) meV. The next-nearest-neighbor inter-chain exchange constant J' is also estimated to 
be 0.05 meV. The observed magnetic excitation spectrum from inelastic neutron scattering is in 
excellent agreement with numerical calculations based on the Mtiller ansatz. 

PACS numbers: 75.10.Pq, 75.25.+Z, 75.50.Ee 
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Progress in the field of quantum magnetism is driven 
by advances in synthesis of prototypical materials. 
One of them is the S=l/2 trimer chain system* 1 * 2 -^ 
One realization of such model was made by Ajiro et 
al. in the material 2Dioxane-3CuCl2^ After follow- 
ing the published sample preparation procedure as de- 
scribed in the Ref.[5], we found this compound to 
be unstable under ambient conditions. Instead, a 
novel compound, 2Dioxane-2H20-CuCl2 (CuDCl) was 
obtained. Here, we describe the crystal structure 
and magnetic property characterization of CuDCl, 
which turns out to be a weakly coupled quasi one- 
dimensional (ID) Heisenberg antiferromagnetic (HAF) 
S=l/2 chain. The physical properties of HAF S=l/2 
chain system have been thoroughly studied in a num- 
ber of materials, such as KCuF 3 6 ^, Cu(C 6 D 5 COO) 2 ■ 
3D 2 (Copper Benzoate^O , BaCu 2 Si 2 ? ii'i2il2ii4il5il£, 
Cu(C4H 4 N 2 )(N03)2 (CuPzN)!i!£. However, unlike 
these compounds, the CI concentration in CuDCl can 
be arbitrarily tuned by Br substitution. As a result, 
the strength of superexchange coupling J changes dra- 
matically with different Br concentrations The new 
compound opens a new research field based on the HAF 
S=l/2 chain system with bond randomness. 

We have determined the crystal structure of CuDCl 
(C 8 H 20 O 6 CuCl2) by X-ray diffraction at T=173 K. The 
compound crystallizes in the monoclinic space group 
C2/c with lattice constants a=17.4298(9) A, 6=7.4770(4) 
A, c=11.8230(6) A, and /3=119.4210(10)°. 

The packing structure can be described as composed 
of CuCi2-2H20 bi-layers in the a-c plane as shown in 
Fig. [Ua). They are well separated by the 1,4-Dioxane 
organic molecules. Hence the possible mediated inter- 
layer interactions existing between molecular units dis- 
placed by (a ± b)/2 should be very weak. Fig. [TJb) 
shows the structure in a single a-c plane. The Cu 2+ 
ions form chains along the crystallographic c axis and 



the distance between nearest neighbor Cu 2+ ions is c/2. 
The chains are separated from each other by the 1,4- 
Dioxane molecules. The Cu 2+ ions have a 4 coordination 
with their neighboring chlorine and oxygen atoms. The 
four short bonds including two Cu-Cl (2.285 A) and two 
Cu-0 bonds (1.959 A) are coplanar, thus defining the 
x-y plane for the 3d x 2_ y 2 copper orbital. Thus the ex- 
change constant coupling the nearest Cu 2+ ions through 
the super-exchange pathway Cu-Cl-H-O-Cu (dash lines) 
should be dominant. 

Crystals of CuDCl were grown by slow evaporation of 
aqueous solutions containing anhydrous Cu(II) chloride 
and excess of 1,4-Dioxane and methanol. Blue single 
crystals with a shape of planar irregular hexagon were 
harvested after a few weeks. 

X-ray diffraction data were collected on a tiny single 
crystal sample using Bruker SMART Apex CCD diffrac- 
tometer (Moka radiation, A=0.71073A) using (j> and ui 
scans. Totally 6122 reflections were collected. The struc- 
ture was solved via direct methods using SHELX97 and 
refined via least-squares using SHELX97. 

DC magnetic susceptibility measurement was per- 
formed on a single crystal sample (m=0.03 g) in the 
temperature range of 1.8 K< T <300 K using a super- 
conducting quantum interference device (SQUID) mag- 
netometer in a DC field of 0.1 Tesla. 

Inelastic neutron scattering (INS) measurements on 
cold neutron triple axis spectrometer TASP^S at the 
SINQ spallation source were performed using one deuter- 
ated single crystal with a mass of m=0.6 g. In this 
Setup, the spectrometer was operated in fully open hori- 
zontal collimation mode with final energy fixed at 5 meV. 
Experiment in Setup 2 was performed on SPINS at the 
NIST Center for Neutron Research (NCNR). The sample 
consists of two co-aligned single crystals with total mass 
0.5 g and a mosaic of 1°. Open horizontal beam diver- 
gences through the triple axis spectrometer were given by 



2 




FIG. 1: (color online) Crystal structure of CuDCl (a) View 
along the c axis showing the well separated bi-layers struc- 
ture, (b) View along the b axis showing how Cu 2+ ions are 
linked into a ID chain. The chain axis is along the c direc- 
tion. Dashed lines show the super-exchange Cu-Cl-H-O-Cu 
pathway. The grey outline is one unit cell. Color coding is as 
follows: Pink:Cu, GreemCl, Blue:C, Red: H, and Yellow:0. 
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FIG. 2: Magnetic susceptibility as a function of temperature 
in a single crystal sample of CuDCl. Inset: enlarged low- 
temperature part of magnetic susceptibility, the solid line is 
a fit to the model as described in the text. 



( 58 Ni-guide)-80'-300'-240' collimations. A horizontally- 
focusing pyrolytic graphite (PG 002) analyzer set to re- 
flect neutrons with final energy Ef — 3.7 meV was used. 
A liquid-nitrogen-cooled Be (Setup 1) or BeO (Setup 2) 
filter was placed after the sample to eliminate higher- 
order beam contamination. In both Setups, sample was 
oriented with its reciprocal (h, 0, 1) plane in the horizon- 
tal plane and the temperature was controlled by the con- 
ventional liquid helium cryostat. Wave- vector transfer is 
indexed as q = ha* + Ic* . 
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FIG. 3: (color online) (a) False-color contour map of mea- 
sured inelastic neutron scattering intensity in CuDCl as a 
function of transferred energy hu) and transferred wavevector 
along the (1,0,1) direction at T=1.6 K, where the background 
is subtracted as described in the text, (b) Instrument reso- 
lution convolved model calculation of two-spinon continuum 
Miiller approximation with J=0.85 meV . Solid lines are the 
predicted lower and upper bounds of the spinon continuum 
with J=0.85 meV. 



The DC magnetic susceptibility x( T)—M/H for a sin- 
gle crystal sample of CuDCl is shown in Fig. [2] The 
inset shows the low-temperature part of x{T)- We ob- 
serve a broad maximum and a rapid decrease at lower 
temperatures, which indicates short range antiferromag- 
netic correlations for low-dimensional systems. The be- 
havior is totally consistent with what is expected for a 
Heisenberg S=l/2 AF chain with a Bonner-Fisher max- 
imum at T max =7.5 K. There is no evidence of 3D long 
range order down to the temperature T=1.8 K from the 
susceptibility data. 

While the susceptibility data is well described by the 
Bonner-Fisher curve, such measurement can not distin- 
guish which direction the chain axis follows. In view of 
the crystal structure, the chain axis is expected to be 
along the crystallographic c axis . To confirm that, we 
measured the magnetic dispersion relations both along 
the chain direction and perpendicular to the chain direc- 
tion in the (h, 0, 1) plane. 

Fig.G3a) shows the false-color map of background sub- 
tracted magnetic scattering intensity for two-spinon con- 
tinuum in CuDCl, at T=1.6 K using Setup 1. This data 
set was obtained by combining eleven constant-q (1,0,/) 
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FIG. 4: Typical scans of background-subtracted neutron scat- 
tering intensity vs transferred energy hu> in CuDCl at T—1.9 
K (a) zone center (h,0,-l) and (b) zone boundary (ft,0,-0.5). 
Solid lines are the fits to the data based on Miiller approxima- 
tion of two-spinon continuum convolved with the instrument 
resolution. 



scans along the chain axis for 0.5< I <1 with step size 
0.05. The background is determined as gaussian fits to 
the raw data which are at least 0.1 meV outside of the 
known bounds of the two-spinon continuum for CuDCl, 
which accounts for the incoherent elastic scattering of the 
sample and the thermal diffuse scattering of the analyzer. 

INS measurements along the chain direction were also 
performed using Setup 2 with better energy resolution. 
Fig. Q] shows the typical background-subtracted neutron 
scattering intensity versus fku. Here the data were col- 
lected by using the focusing analyzer mode, where the 
value of h was chosen so that the scattered wave-vector 
kf is parallel to the chain axis c to achieve the good 
wave- vector resolution along the chain. Fig. 01a) at zone 
center q=(/i,0,-l) shows the scattering intensity consis- 
tently decreasing, which indicates the divergence at the 
lower bound in the continuum. Fig. EJb) at zone bound- 
ary q=(/i,0,-0.5) shows a prominent peak at ~1.4 meV, 
which is corresponding to the value of ^ and J is the 
nearest neighbor intra-chain exchange constant. 

Furthermore, we measured the dispersion perpendic- 
ular to the chain axis using Setup 2. Representative 
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Setup 2 
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FIG. 5: Typical scans of background-subtracted neutron scat- 
tering intensity vs hoj in CuDCl at T—1.9 K showing the weak 
dispersion of magnetic excitation with h. Solid lines are the 
fits to the data as described in the text. 



background-subtracted energy scans up to 1 meV with 
wave- vectors perpendicular to the chain axis c are shown 
in Fig. [5] To determine the background, constant-q scans 
of (/i,0,-0.5) and (/i,0,-1.5) at the same configurations 
were measured for each h, where no magnetic excitation 
is expected in this energy transfer range (Hu> <1 meV) 
and the background is taken account as a fit to Gaussian 
profile. The apparent dispersion along h suggests the 
existence of next-nearest-neighbor exchange constant J' 
along the a axis. 

The solid line in Fig. [5] is a fit of measured DC suscep- 
tibility data to 



X = Xdia + f X Ximp + (1 - /) X X 



(1) 



where \dia is the temperature-independent diamagnetic 
contribution of the sample, Ximp is the contribution from 
the impurity, x-m is the spin susceptibility for ID S=l/2 
AFM chain and / is the fraction of free spin due to the 
impurity. Ximp simply obeys the Curie law. For Xm, a 
polynomial expression^ which can reproduce the results 
of Bonner and Fish22 very well was used. 

The fitting procedure yields Xdia=-2A(1) x 10 -4 
emu/Cu mol, /=0.55(2)%, and J=0.94(3) meV. 

Since the numerical difference for the dynamic spin- 
correlation function <S(q, to) between the exact two- 
spinon solution^-s^ 4 . and Miiller approximation^ is small, 
the data in Fig. [3)Ja) were globally fitted to the Miiller 
approximation of the ID HAF S=l/2 chain convolved 
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with instrumental resolution: 

s( )= 1 A 

[q,UJ> 2^(^)2 _ e?(q) 
X@(fuv - ei(q))6(e 2 (q) - hu), 



(2) 



where A is an overall factor, Q(x) is the step func- 
tion, ei(q)=^| sin(7rZ)| and e 2 (q)=7rJ| sin ^-\ arc the 
lower and upper bounds of the two-spinon continuum. 
Fig. [2Kb) shows the best fit of model simulation with 
J=0.85(4) meV. The intensity is high at the lower bound 
due to the inverse square root divergence of the two- 
spinon density of states at that boundary. The solid lines 
are the lower and upper two-spinon continuum bound- 
aries with J=0.85 meV. 

To estimate the inter-chain coupling J', we need to 
consider the case of weakly coupled AFM S=l/2 chains. 
Although the spin dynamics was well studied in the or- 
dered phase by using the mean-field (MF) approxima- 
tion and random-phase-approximation (RPA) 26 ' 27 , un- 
fortunately so far there is no such theory to predict the 
spin dynamics including the transverse dispersion above 
Tjy. Nevertheless, we naively treat the transverse disper- 
sion as a simple sinusoidal function. The lower bound 
of the continuum is modified to also include the disper- 
sion along the transverse direction in the same way as 
the transverse dispersion term is included in spin wave 
below TV: 



7T 2 J 2 

e 2 (q) = — — sin 2 (7rZ) + A 2 sin 2 (7r/i), 



(3) 



where A is the bandwidth of the transverse dispersion. 
The expression of upper bound of the continuum is un- 



changed. The solid lines in Fig. [4] and Fig. [5] are the fits 
to the instrument convolution of the model with fixed 
J=0.85 meV and the only adjustable parameter is the 
overall factor A. It agrees very well with the data and 
the fitting procedure gives A=0.29(l) meV. We further 
assume that the bandwidth A is related to J' as the same 
way in the transverse dispersion below Tjy^, J' can be 
roughly estimated to be ~A/6=0.05 meV. Using the val- 
ues of J and J', the Eq. (13) in the Ref.[26] gives the 
estimation of Tm ~1.3 K. 

In summary, we have carried out x-ray crystallogra- 
phy, magnetic susceptibility and inelastic neutron scat- 
tering measurements to understand what is the right spin 
Hamiltonian of single-crystal samples in a newly syn- 
thesized compound CuDCl. In view of crystal struc- 
ture, it consists of ID S=l/2 HAF chains with the 
chain axis along the c axis. Bulk measurement and in- 
elastic neutron scattering data are also consistent with 
the weakly coupled chain model calculation. Our fu- 
ture work will focus on the mixed magnetic system 
2Dioxane-2H20-Cu(Cli_ 2: Br 2: )2, which could be an ideal 
realization of S=l/2 chain with bond randomness that 
to date has only been theoretically studied] 28 ' 29 
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